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ABSTRACT 

A  computer  model  is  presented  for  determining  the  daily  allo- 
cation of  airstrikes  for  interdicting  a  lines-of-communication  (LOC) 
network  assuming  an  exponential  damage  function.   The  allocation  of 
sorties  to  arcs  in  the  network  is  based  on  the  assumptions  that  (1) 
flow  of  supplies  is  restricted  by  network  capacity  and  (2)  parameter 
values  for  each  arc  are  independent.   Upper  and  lower  flow  bounds, 
cost  per  unit  flow,  maximum  repair  time  and  cost,  and  the  vulner- 
ability parameter  for  each  arc  are  required  as  input  data.   The 
model  selects  that  arc  to  strike  which  maximizes  the  repair  cost 
plus  the  product  of  the  increase  in  minimum  cost  circulation  flow 
and  repair  time.   The  procedure  is  programmed  in  daily  cycles, 
allowing  repair  m  interdicted  arcs.  A  sample  problem  and  all  docu- 
mentation necessary  for  duplicating  the  computer  program  are  given. 
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I.   INTRODUCTION 

A.    GENERAL 

The  ability  of  a  military  force  to  conduct  sustained  combat 
operations  any  distance  from  its  base  of  supply  may  be  limited  by 
the  land  lines-of-communication  (LOC)  network  over  which  its 
supplies  must  flow.   This  is  particularly  true  if  air  and  sea 
resupply  are  not  feasible,  or  the  opposing  force  has  superiority 
over  these  alternatives.   The  amount  and  type  of  material  required 
to  support  a  military  force  are  determined  by  the  nature  of  combat 
operations  and  the  intensity  of  combat.   A  large  conventional  force 
would  require  more  pounds-per-man-day  than  a  small  guerilla  force 
capable  of  obtaining  a  portion  of  its  requirements  in  the  area  in 
which  it  operates.   Undoubtedly,  however,  a  requirement  will  exist 
for  supplies  in  excess  of  what  can  be  obtained  locally  or  what  can 
be  transported  organically  by  the  combat  force. 

The  land  LOC  network  then  becomes  an  important  military  target. 
Flow  of  military  material  over  the  network  can  be  reduced  by  any  of 
three  alternatives:   (1)  attacking  the  transport  capability,  i.e., 
vehicles  or  rolling  stock,  (2)  attacking  manufacturing  and  supply 
centers,  thereby  reducing  the  amount  of  materiel  available,  or  (3) 
attacking  the  road/rail  LOC  network  over  which  the  materiel  must  flow. 
Certainly,  a  combination  of  these  alternatives  would  be  most  desirable, 
However,  political  considerations  might  preclude  attacking  manufactur- 
ing and  supply  centers;   and  terrain,  weather,  and  lack  of  intell- 
igence information  might  limit  the  effectiveness  of  attempting 


to  attack  convoys  moving  over  the  network. 

Such  a  scenario  is  particularly  applicable  to  the  present  war 
in  Southeast  Asia  and  is  a  general  characteristic  of  limited  warfare. 
The  North  Vietnamese  neither  possess  a  sophisticated  air  resupply 
capability,  nor  can  they  openly  move  supplies  by  sea  to  support 
their  forces  operating  in  South  Vietnam.   They  do,  however,  enjoy 
immunity  from  attack  on  all  their  facilities  north  of  the  17th 
Parallel.   Therefore,  the  supply  of  materiel  to  their  forces  oper- 
ating in  the  South  is  primarily  limited  by  the  transportation  system 
over  which  this  materiel  must  flow. 

The  objective  of  an  interdiction  campaign  by  airstrike  must 
then  be  to  systematically  reduce  the  capacity  of  road  or  rail 
segments  which  will  ultimately  restrict  the  amount  of  materiel  that 
can  flow  through  the  network. 

B.    AVAILABLE  LITERATURE 

Several  authors  have  published  works  on  the  problem  of  inter- 
dicting a  lines-of-communication  network.   Durbin  [3]   has  construct- 
ed a  model  for  determining  maximum  cargo  flow  as  a  function  of 
available  transport  vehicles.   His  method  sequentially  selects  and 
interdicts  the  most  vital  segment  in  the  network  until  flow  is  re- 
duced to  zero  or  a  specific  number  of  segments  have  been  destroyed. 
The  model  operates  in  daily  cycles  and  allows  repair  on  interdicted 
segments. 

Wollmer  [11,  14,  15]  has  developed  a  procedure  to  determine 
the  n  most  vital  segments  in  a  network  and  to  evaluate  the  effect 


of  removing  the  most  vital  segment.   In  the  latter  work,  road 
capacity  as  a  function  of  the  number  of  strikes  is  reduced  (1) 
by  a  known  quantity  and  then  (2)  by  a  quantity  that  is  a  random 
variable. 

Wollmer  [12,  13]  has  also  developed  algorithms  and  a  model 
for  evaluating  a  targeting  strategy  whose  aim  is  the  greatest  re- 
duction in  network  flow.   The  model  is  programmed  in  daily  cycles 
with  the  user  specifying  the  number  of  days  and  strikes. 

Mustin  [9]  uses  the  topological  dual  to  represent  a  lines-of- 
communication  network  and  assumes  the  decrease  in  capacity  of  a 
road  or  rail  segment  has  a  fixed  deterministic  rate  which  is  linear 
between  its  upper  and  lower  bounds. 

Nugent  [10]  also  uses  the  topological  dual  for  network  rep- 
resentation, but  assumes  th/:  reduction  in  capacity  of  a  road  or 
rail  segment  is  exponential  with  decreasing  marginal  returns .   His 
model  is  constrained  by  the  number  of  available  sorties  and  does 
not  allow  for  the  dynamics  of  the  daily  allocation  problem. 

C.    ALLOCATION  OF  EFFORT 

Many  techniques  have  been  devised  to  aid  in  solving  the 
problem  of  allocating  limited  resources.   Hancock  [8]  provides 
a  reference  on  classical  optimization  techniques;   Danskin  [1,  2] 
considers  the  theory  of  maxima  and  minima  and  its  application  to 
weapons  systems  allocation;   and  Koopman  [6,  7]  considers  a  simi- 
lar problem  in  the  allocation  of  search  and  screening  effort. 
The  problem  facing  the  strike  planner  is  to  determine  what  road 


or  rail  segments  in  the  network  to  attack  and  how  many  strikes  to 
allocate  to  each  selected  segment.   Numerous  considerations  must 
be  weighed  in  order  to  make  this  determination:  (1)  the  number  of 
sorties  available  on  a  given  day,  (2)  the  probability  that  a  strike 
against  a  particular  segment  is  successful,  (3)  the  expected  cost 
to  the  enemy  in  time  and  effort  to  repair  an  interdicted  segment, 
(4)  the  expected  cost  of  aircraft  lost  in  striking  a  defended  seg- 
ment, and  (5)  the  expected  reduction  in  capacity  if  a  strike  is 
successful. 

The  obvious  procedure  in  making  the  allocation  of  W  sorties 
is  to  attack  those  segments  (at  a  level  of  n. .  each)  that  maximize 
the  ratio  of  benefits  to  cost.   As  will  be  shown  later,  Gibbs' 
Lemma  [2]  provides  the  mathematical  basis  for  determining  this 
allocation. 
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II.   OBJECTIVE  AND  SCOPE 

The  objective  of  this  paper  is  to  present  a  model  for  determin- 
ing the  daily  allocation  of  airstrikes  for  interdicting  a  lines-of- 
communication  network  assuming  an  exponential  damage  function.   The 
model  seeks  to  allocate  a  maximum  of  W  (specified)  sorties  per  day 
to  those  arcs  (segments)  which  yield  the  greatest  reduction  in  flow 
capacity  and  the  greatest  time  and  cost  to  repair.   Two  options  are 
provided  that  will  allow  the  strike  planner  to  (1)  specify  the 
number  of  strikes  per  day  and  the  number  of  days  in  the  interdiction 
campaign,  or  (2)  select  a  percentage  of  the  maximum  possible  flow 
to  be  reduced.  Using  the  second  option,  the  maximum  available 
number  of  sorties  per  day  and  the  number  of  campaign  days  are 
specified. 

Certain  information  is  assumed  to  be  available  to  the  strike 
planner  for  each  arc  in  the  LOC  network: 

(1)  Upper  flow  bound  -  maximum  capacity  in  units  per  day. 

(2)  Lower  flow  bound  -  minimum  capacity  in  units  per  day  below 
which  the  flow  cannot  be  reduced. 

(3)  Maximum  repair  time  -  the  time  required  to  repair  the 
maximum  possible  damage  on  an  arc. 

(4)  Maximum  repair  cost  -  the  cost  to  repair  the  maximum 
possible  damage  on  an  arc. 

(5)  Vulnerability  parameter  -  the  probability  that  a  sortie  is 
successful.   The  strike  planner  can  manipulate  this  para- 
meter by  providing  his  own  definition  of  "successful." 

The  vulnerability  of  any  arc,  V..,  is  defined  as  the  upper  flow 
bound  minus  the  lower  flow  bound,  i.e.,  U  -L   .   A  necessary  assump- 
tion is  that  the  vulnerability,  repair  cost,  and  repair  time  are  in- 
dependent of  the  same  parameters  for  any  other  arc. 
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By  assuming  an  exponential  damage  function,  the  reduction  in 
capacity  d,.,  for  any  arc  (i,j)  caused  by  n. .  strikes  is: 

dij  =  vijri "  ^p^u  n±pi  <x> 

Arc  repair  time  and  cost  are  linear  functions  of  the  amount 
of  damage  caused  by  n. .  strikes: 

t«-«yvi«  (2) 

hi  •  (WKij  (3) 

The  amount  of  arc  repair  that  can  be  completed  in  any  one  day 
is  a  function  of  the  computed  cumulative  arc  repair  time  and  the 
total  damage  on  the  arc: 

r. .  =  (1/CRT. .)  D..  (4) 

The  value  of  allocating  n. .  strikes  per  day  to  arc  (i,j) 


is  determined  by: 


Value. .  =  (F'  -  F)  t. .  +  k. .  (5) 

3-3  ij    iJ 


The  model  daily  selects  and  interdicts  those  arcs  that 
maximize  (5) ,  performs  repair  as  indicated  by  (4) ,  and  steps  to 
the  next  day.   The  procedure  has  been  programmed  in  FORTRAN  IV 
for  use  on  the  IBM  360/67  computer. 
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III.   NETWORK  REPRESENTATION  OF  LINES-OF-COMMUNICATION 

A.    NETWORK  STRUCTURE 

A  transportation  system,  whether  single  or  multimode,  can  be 
represented  by  a  connected  graph  of  nodes  and  directed  arcs.   For 
example,  a  highway  system  would  consist  of  roads  (arcs)  and  the 
intersection  of  roads  (nodes) .   If  a  road  is  one-way,  the  beginning 
node  of  an  arc  is  the  point  where  traffic  enters  the  segment  and  the 
ending  node  is  the  point  where  traffic  leaves  the  segment.   If  a 
road  is  two-way,  it  is  represented  by  two  arcs,  one  in  each  direction. 

The  LOC  user  attempts  to  send  units  of  flow,  expressible  in 
vehicles  per  day,  tons  per  day,  etc.,  from  the  source  to  the  sink. 
The  source  of  flow  could  be  a  rear-area  supply  center  and  the  sink 
could  be  the  battlefront  or  a  forward-area  supply  center.   If  the 
user  is  constrained  by  the  capacity  of  the  transportation  system, 
it  is  assumed  that  he  will  attempt  to  maximize  the  total  flow  of 
cargo  through  the  network  at  a  minimum  cost. 

Consider  the  example  LOC  networ  shown  in  Fig.  1.   Suppose  the 

user  wishes  to  send  flow  from  node  AA  to  node  ZZ.   Node  AA  is  the 

source  and  node  ZZ  is  the  sink.   The  route  (AA,A,K,L,ZZ)  is  one  of 

many  possible  routes  over  which  the  flow  can  be  sent.   Associated 

with  each  arc  are  certain  parameters:   lower  bound  on  flow,  upper 

bound  on  flow  (flow  capacity),  and  the  cost  per  unit  flow.   For  arc 

(K,L)  the  following  parameters  are  given: 

Lower  bound  (L  )   0 

Upper  bound  (U  )   400 

Flow  cost  (C  )     7 
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These  parameters  imply  that  between  0  and  400  units  of  flow  per  time 
unit  can  be  sent  over  (K,L)  at  a  cost  of  7  per  unit.   This  cost  can 
be  expressed  in  terms  of  dollars  per  ton  of  cargo,  dollars  per 
vehicle,  manpower  units,  or  any  other  definable  cost.   If  the  arc 
is  used  to  its  capacity,  i.e.,  M..  =  U..,  it  would  cost  2800  per  time 
unit  to  send  flow  over  this  arc.   The  total  cost  per  time  unit  of 
sending  400  units  of  flow  from  node  AA  to  node  ZZ  over  route 
(AA,A,K,L,ZZ)  would  be: 

I  M..C..  =  400(0)  +  400(6)  +  400(7)  +  400(0) 

all  (i,j)   1J  13 

in 
the  route        =  4200 

It  is  important  to  note  that  the  flow  out  of  a  node  equals  the 
flow  into  a  node.   This  conservation  of  flow  implies  there  can  be  no 
storage  at  nodes.   Therefore,  for  most  general  networks,  if  X  units 
of  flow  are  sent  from  the  source  to  the  sink,  then  any  flow  pattern 
that  accomplishes  this  will  have  net  flows  of  X  units  out  of  the 
source,  X  units  into  the  sink,  and  zero  on  all  other  nodes. 

In  actual  LOC  networks,  there  may  be  more  than  one  source  and/or 
sink.   This  can  be  handled,  provided  any  source  can  supply  any  sink, 
by  adding  artificial  nodes  and  arcs.   This  allows  the  user  to  con- 
vert the  multi-source,  multi-sink  problem  to  a  single-source,  single- 
sink  problem.   The  artificial  arcs  are  used  to  connect  the  true 
sources  (sinks)  to  an  artificial  super  source  (sink) .   Note  in 
Figure  1  that  arcs  (AA,A) ,  (AA,B),  (AA,C),  and  (AA,D)  can  be  con- 
sidered artificial  arcs  connecting  the  true  sources  A,  B,  C,  and  D 
to  the  single  artificial  super  source  AA.   Similarly,  arcs  (L,ZZ), 


15 


(H,ZZ),  (I,ZZ),  and  (J,ZZ)  can  be  considered  artificial  arcs  connect- 
ing the  true  sinks  L,  H,  I,  and  J  to  the  artificial  super  sink  ZZ. 
The  flow  on  an  artificial  arc,  e.g.,  (AA,A),  can  be  considered  the 
rate  that  cargo  leaves  the  true  source  A.   The  capacity  for  arti- 
ficial arcs  is  usually  infinite,  but  could  represent  the  maximum 
supply  rate  at  a  facility  such  as  a  supply  center  or  warehouse. 
The  lower  flow  bound  and  cost  on  artificial  arcs  are  normally  zero. 
Artificial  arcs  are  not  considered  vulnerable  to  attack. 

Actual  flow  on  an  arc  is  designated  M. .  if  the  flow  is  from 
node  i  to  node  i  and  M..  if  the  flow  is  in  the  reverse  direction. 
Normally,  flow  is  assumed  to  be  from  the  source  to  the  sink  but 
can  be  in  either  direction  on  intermediate  two-way  arcs. 

U . . ,  U..,  and  L. .  represent  the  upper  flow  bound,  the  flow 

capacity  at  any  time,  and  the  lower  flow  bound,  respectively,  for 

any  arc  (i,j).   The  following  relationship  must  hold  for  these 

parameters: 

0  <  L. .  <  M. .  <  U. .  <  U#J 
-  ij  -  ij  -  ij  -  lj 

The  maximum  possible  flow  through  the  network  can  be  determined 

by  using  the  maximal-flow  minimal-cut  theorem  14].   It  states  that 

the  maximum  flow  in  a  network  is  equal  to  the  minimum  value  of  all 

cut-sets,  where  a  cut-set  is  defined  as  a  set  of  arcs  that  separate 

the  source  from  the  sink.   The  value  of  a  cut-set  is  defined  as  the 

sum  of  the  capacities  of  all  arcs  in  the  cut-set.   In  the  example 

network,  Figure  2,  arcs  (A,K),  (K,E),  (E,H),  (E,F),  (C,F),  (C,G), 

and  (D,G)  form  the  minimum  cut-set  with  a  value  of  5000  units. 
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B.    CIRCULATION  FLOWS  IN  NETWORKS 

As  previously  stated,  if  the  LOC  user  is  constrained  by  the 
capacity  of  the  network,  it  is  assumed  he  wishes  to  maximize  the 
total  flow  through  the  network  at  a  minimum  cost.   Such  a  problem 
is  efficiently  handled  by  the  "Out-of-Kilter"  algorithm  developed 
by  Fulkerson  [4].   The  network  is  restructured  by  adding  a  "universal" 
arc  from  the  super  sink  to  the  super  source  as  in  Fig.  3.   This 
artificial  arc  is  assigned  an  upper  flow  bound  greater  than  or  equal 
to  the  maximum  possible  f low  through  the  network  (usually  infinite) , 
a  lower  bound  of  zero,  and  a  unit  flow  cost  less  than  the  flow  cost 
on  any  arc  (usually  minus  infinity).   Hence,  the  problem  of  maxi- 
mizing source-sink  flow  is  converted  to  one  of  finding  a  minimum 
cost  circulation  flow  by  maximizing  flow  on  the  universal  arc. 


18 


19 


IV.   MODEL  OPERATION 


A.    MATHEMATICAL  PRELIMINARIES 


Assume:   (1)  The  set  X*  =  (x*  x*  ,...,  x*) 

1   L  n 


maximizes    E  f  .  (x .  ) 


n 

I 

i=l 


n 
sub i .to      £  x .  =  X   ;x.    0 
i-1  X 

(2)  The  f.(x.)  are  continuous  and  diff erentiable 

(3)  f\'  (x.)  <  0 

Note:   These  are  necessary  and  sufficient  conditions 
for  a  global  maximum  since  the  objective  func- 
tion is  concave  and  the  constraints  are  convex. 

Then  there  exists  a  scalar  s  >  0  such  that: 

(1)  f!(x.)  =  s  for  x*  >  0 

(2)  f1(x.)  <  s  for  x*  =  0 

11  l 

Note:   The  Lagrange  Multiplier,  s,  can  be  interpreted 
as  the  implicit  marginal  price  of  allocating  an 
additional  unit  of  resource. 

The  above  is  known  as  Gibbs1  Lemma  and  is  discussed  in  more  detail  by 

Danskin  [2].   It  can  be  derived  from  the  Kuhn-Tucker  conditions  or  by 

elementary  arguments  based  on  the  special  structure  of  the  problem. 

There  are  two  important  implications  of  this  lemma: 

(1)  Marginal  return  must  be  above  a  threshold  value,  s,  for  a 
resource  to  be  allocated 

(2)  For  alternatives  which  receive  allocation,  f.(x.)  =  f.(x.), 

liii 
i.e.,  at  the  optimum  allocation  point,  the 

ratio  (marginal  return/marginal  cost)   is  constant. 
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B.    MODEL  DESCRIPTION 

As  stated  earlier,  there  are  numerous  considerations  to  be 
weighed  in  allocating  strikes  to  arcs  in  the  LOC  network.   From  a 
mathematical  viewpoint,  it  is  assumed  that  the  effectiveness  of  the 
network  to  its  user  may  be  measured  in  terms  of  how  to  obtain  maxi- 
mum flow  at  minimum  cost.   Therefore,  the  criterion  to  be  used  to 
allocate  sorties  is  to  minimize  LOC  effectiveness  over  time,  taking 
into  account  the  cost  to  repair  damage  on  an  arc.   This  is  accomplish- 
ed by  allocating  strikes  to  the  set  of  arcs  that  maximize  the  costs 
to  the  user. 

The  effect  of  striking  an  arc  is  to  reduce  its  capacity  and 
increase  the  cost  of  a  minimum  cost  circulation  flow  for  a  specified 
period  of  time.   Therefore,  in  allocating  hi   strikes,  the  selected 
arcs  are  the  ones  that  maximize  repair  cost  plus  the  product  of  the 
increase  in  minimum  cost  circulation  flow  and  repair  time.   Specific- 
ally, if  F  is  the  cost  of  a  minimum  cost  circulation  flow  before  arc 
(i,j)  is  struck,  F  is  the  cost  of  a  minimum  cost  circulation  flow 

after  n. .  strikes  on  arc  (i,i),  and  t..  and  k. .  are  the  repair  time 
xj  xj       1J 

and  repair  cost,  respectively,  for  damage  caused  by  n. .  strikes, 

then  the  arcs  selected  for  attack  will  be  the  ones  that  maximize: 

£     (F*  -  F)  t..  +  k,.  " 
all   (i,j)  lj    1J 

s.t.    (1)     I     n   =  N   ;  n. .  >  0 
all  (i,j)  iJ         1J 

where:  F  =    Z     M  C.. 

all  (i,j)  i:>  lJ 
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F  -    I  M..C..   -  d..(C..  +  Cu)   ;  C  >  0 

all  (i,j)  1J  lj     ±J   ij  1J 

CU  «  0 

Note:   C  is  the  unit  flow  cost  on  the  universal 
arc, 

dij=VijIl-exp(-VPijnij)] 

t..  =  (d../V..)  T. . 

k..  =  (d../V..)  K.. 

The  objective  function  then  becomes  to  maximize: 

-.E,,  ..-T.  .V.  .(C.  .  +  CU)  [1  -  exp(-VP.  .  n.  ,)]2 
all  (i,j)   iJ  iJ   i]  ij  ij 

+  K. .11  -exp(-VP. .  n. .)] 

This  can  be  shown  to  be  a  continuous  concave  function  of  n. .  whose 
second  derivative  is  negative.   The  constraint  is  convex;   therefore, 
the  above  objective  function  becomes  the  f.(x.)  of  Gibbs'  Lemma. 

To  determine  the  allocation  of  W  airstrikes,  it  is  necessary  to 
compute  the  value  of  the  Lagrange  Multiplier,  s  .   By  Gibbs  Lemma, 
it  is  known  that  there  exists  an  s*  >  0  such  that: 

n#J  =  0  for  f  !(X.)    <  s* 

n.  .  >  0  for  f '  (xJ)    >  s* 
ij  1  i 

An  iterative  method  for  this  computation  is  developed  by  Nugent  [10] . 

The  above  procedure  is  used  in  Option  (1)  where  the  number  of 

strikes  is  specified.   For  Option  (2),  a  similar  procedure  is 

developed  with  the  addition  of  another  constraint.   Recall  that 

Option  (2)  requires  a  percentage  by  which  the  initial  flow  is  to 
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be  reduced.   The  additional  constraint  becomes: 

(2)   I  M   >  (1  -  P)m! 

all  (i,j)  1J  -         1J 

i 

where:  M. .  is  the  initial  network  flow. 

p   is  the  percentage  by  which  the 
initial  flow  is  to  be  reduced. 

Obviously,  for  the  allocation  of  a  maximum  of  N   sorties,  either 

constraint  (1)  or  (2)  will  be  active  and  the  value  of  the  Lagrange 

Multiplier,  s  ,  can  be  computed  using  Nugent  s  method. 

If  the  vulnerability  parameter  of  arc  (i,j)  takes  cognizance 

of  the  probability  of  losing  an  aircraft  on  any  mission  against 

that  arc,  then  the  model  seeks  to  interdict  that  set  of  arcs  where 

the  ratio  (marginal  return/marginal  cost).,  is  constant. 

C.    COMPUTER  SOLUTION  PROCEDURE  AND  OUTPUT 

The  computer  solution  procedure  operates  in  daily  cycles,  the 
strike  planner  specifying  by  Option  (1)  or  (2)  the  specified  or 
maximum  number  of  strikes  per  day,  respectively,  and  the  number  of 
days  in  the  interdiction  campaign.   The  computer  solution  differs 
from  the  optimal  model  described  in  B.,  in  that  strikes  are  allocated 
one  at  a  time.   Each  time,  the  arc  selected  to  be  struck  is  the  one 
which  maximizes: 

(f'  -  F)  t. .  +  k..  . 

i 
The  values  F  ,  F,  t..,  and  k. .  are  computed  for  one  strike  instead 

of  for  n. .  strikes.   This  is  a  marginal  allocation  procedure,  the 
13 

solution  having  the  property  that  if  an  arc  is  selected  for  inter- 
diction when  W  strikes  are  allocated,  it  will  also  be  selected  for 
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interdiction  when  more  than  W  strikes  are  allocated.   Since  optimal 
allocation  of  strikes  would  not  necessarily  have  this  property,  the 
selected  arcs  are  optimal  for  one  strike,  but  only  approximate  optim- 
ality  for  multiple  strikes. 

After  each  strike,  the  following  information  is  output: 


(1 

(2 
(3 
(4 
(5 
(6 
*(7 
(8 


Arc  selected  for  strike 

Arc  capacity  before  and  after  strike 

Actual  arc  flow  before  and  after  strike 

Day  arc  repair  is  to  be  completed 

Cumulative  repair  cost 

Network  flow  after  strike 

Network  flow  cost  after  strike 

Number  of  strikes  against  that  arc  for  that  day 


*Note:   Since  the  unit  flow  cost  on  the  universal  arc  is  an 
artificial  cost  used  to  convert  a  particular  type  of  source- 
sink  flow  problem  to  one  of  circulation  flow,  it  is  not  con- 
sidered in  the  printed  output  calculation. 

The  end  of  a  campaign  day  is  determined  by  the  option  used. 
Using  Option  (1) ,  the  day  terminates  when  W  sorties  have  been 
allocated.  Using  Option  (2),  the  day  terminates  when  either  W 
strikes  have  been  allocated  or  the  total  network  flow  has  been  re- 
duced to  a  specified  percentage  of  the  initial  flow.   The  following 
information  is  output  at  the  end  of  each  day: 

(1)  Arcs  selected  for  strike 

(2)  Number  of  strikes  allocated  against  arc  (i,j)  for 
that  day 

(3)  Number  of  strikes  allocated  against  arc  (i,j)  for 
the  campaign 

(4)  Total  number  of  strikes  allocated  that  day 

(5)  Total  number  of  strikes  allocated  for  the  campaign 

At  the  start  of  each  day,  all  arc  repair  scheduled  for  that 
day  has  been  completed.   The  resulting  flow  pattern  is  computed  and 
arc  flow,  maximum  network  flow,  and  network  flow  cost  are  output. 
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V.    SAMPLE  PROBLEM 

The  network  shown  in  Fig.  4  is  used  as  a  sample  interdiction 
problem.   It  is  the  same  network  used  previously  but  with  the 
addition  of  maximum  repair  time  and  maximum  repair  cost  parameters 
for  each  arc.   Figure  5  depicts  the  same  network,  but  node  positions 
in  relation  to  each  other  have  been  changed  to  indicate  a  more 
realistic  lines-of-communication  network.  Artificial  arcs  from 
the  super  source  are  used  to  indicate  that  flow  in  the  enemy  sanctu- 
ary area  is  invulnerable  to  attack.   Nodes  A,  B,  C,  and  D  are  points 
on  supply  routes  where  arc  capacity  is  subject  to  interdiction.   The 
enemy  objective  is  to  move  supplies  from  his  sanctuary  area  to  for- 
ward supply  centers  located  at  nodes  L,  H,  I,  and  J.   Nodes  K,  E,  F, 
and  G  are  route  intersections.   Figure  5  also  gives  the  actual 
direction  and  magnitude  of  flow  before  the  interdiction  campaign 
begins . 

The  objective  of  the  strike  planner  is  to  reduce  the  initial 
maximum  network  flow  by  60  percent  and  thereafter,  maintain  this 
reduced  flow  for  a  campaign  of  5  days  duration.   He  has  a  maximum 
of  50  aircraft  available  each  day  to  fly  interdiction  missions. 

The  computer  solution  to  this  problem  is  given  in  Appendix  D. 
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ENEMY  TERRITORY 
(sanctuary  area) 


ENEMY  TERRITORY 
(free  fire  area) 


Figure  5.   Network  for  Sample  Problem  Giving  Direction 

and  Magnitude  of  Initial  Flow 
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VI.   SUMMARY  AND  EXTENSIONS 

A  computer  model  has  been  presented  for  determining  the  daily 
allocation  of  airstrikes  for  interdicting  a  land  lines-of- 
communication  (LOC)  network,  assuming  an  exponential  damage  function. 
The  allocation  of  sorties  to  arcs  in  the  network  is  based  on  the 
assumptions  that  (1)  flow  of  supplies  is  limited  by  network  capacity, 
and  (2)  parameter  values  for  each  arc  are  independent.   The  model 
selects  that  arc  to  strike  which  maximizes  the  repair  cost  plus  the 
product  of  the  increase  in  minimum  cost  circulation  flow  and  repair 
time.   The  procedure  operates  in  daily  cycles,  performing  repair  on 
arcs  that  have  been  struck. 

Two  options  are  provided  that  allow  the  strike  planner  to 
(1)  specify  the  number  of  strikes  per  day  and  the  number  of  days 
in  the  campaign,  or  (2)  prescribe  a  percentage  by  which  the  initial 
flow  is  to  be  reduced,  the  maximum  number  of  strikes  per  day,  and 
the  number  of  campaign  days.   The  upper  and  lower  flow  bounds, 
cost  per  unit  flow,  maximum  repair  time  and  cost,  and  the  vulner- 
ability parameter  for  each  arc  in  the  network  are  required  as  in- 
put data. 

The  model  has  been  programmed  in  FORTRAN  IV  for  use  on  the 
IBM  360/67  computer.   Detailed  information  for  each  day  in  the 
campaign  is  listed  in  the  computer  output  to  include  arcs  selected 
for  interdiction,  network  flow  and  network  flow  cost  at  the  end 
of  each  day,  and  the  number  of  sorties  allocated  to  selected  arcs. 
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There  are  almost  an  endless  number  of  extensions  that  can  be 

considered  in  the  allocation  of  airstrikes  to  a  LOC  network.   In 

the  present  model,  it  is  assumed  that  only  one  type  of  aircraft  is 

available  to  attack  the  network,  i.e.,  one  weapons  system  against 

y   targets.  A  more  practical  and  interesting  approach  is  to  assume 

that  more  than  one  type  of  aircraft  is  available  for  interdiction, 

each  with  a  different  kill  probability  against  different  types  of 

targets,  i.e.,  X  weapons  systems  against  y   targets.   The  damage  on 

an  arc  caused  by  n. .  strikes  then  becomes: 

ij 

diJ-ViJ  I1-j1erp(-VP«,1cniJ,k)] 

Another  modification  would  be  to  vary  the  damage  function  it- 
self.  The  exponential  was  used  in  this  model  because  it  exhibits 
decreasing  marginal  returns  to  scale  and  simplifies  computational 
procedure;   however,  other  such  functions  could  be  utilized.   The 
exponential  is  a  continuous  approximation,  for  n. .  large  and  VP.. 


small,  of  the  function: 


v±jn  -  a  -  yp1;J)nijj 


A  contrast  could  be  made  of  the  damage  caused  by  both  distributions 
with  actual  damage  data.   Holliday  [5]  has  done  research  in  the 
field  of  damage  assessment  for  the  Advanced  Research  Projects 
Agency  and  his  work  might  be  used  for  such  a  contrast. 
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Card 

Column 

Form, 

no. 

no. 

1 

1-80 

A 

2 

1-5 

15 

6-10 


11-15 


16-20 


22-26 


15 


15 


15 


F4.2 


Arc 

1-4 

A4 

Data 

Cards 

6-7 

A2 

9-10 

A2 

12-20 

19 

22-30 

19 

32-40 
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APPENDIX  A 
COMPUTER  PROGRAM  INPUT  DATA 
tt  Description 

Program  title.  (PRONAM) 


Maximum  number  of  strikes  that  can  be 
allocated  each  day.  (N) 

The  number  of  days  in  the  campaign. 
(NUMDAY) 

The  number  of  invulnerable  or  artificial 
arcs  in  the  network.   Include  all  arcs 
considered  invulnerable  to  airs tr ike, 
artificial  arcs  connecting  the  dummy 
source  to  the  true  source (s),  and  arti- 
ficial arcs  connecting  the  true  sink(s) 
to  the  dummy  sink.  (INVL) 

If  equal  to  1,  specifies  one-way  flow  on 
some  arcs  in  the  network.   If  equal  to 
2,  specifies  two-way  flow  on  all  arcs  in 
the  network.   When  arcs  are  two-way,  a 
strike  against  one  implies  a  strike 
against  the  other.   (X2WAY) 

Specifies  the  percentage  of  flow  to  be 
interdicted.   If  greater  than  1.0,  a 
fixed  number  of  strikes  per  day  (N)  will 
be  targeted.  (REDU) 

Name  of  arc.   If  two-way  flow,  must  be 
the  same  for  both  arcs.  (ARNAM) 

Beginning  node.   (FROM) 

Ending  node.   (TO) 

Lower  bound  on  arc  flow.   (LBND) 

Arc  capacity  (upper  bound  on  arc  flow) • 
(MCTY) 

Cost  of  one  unit  of  flow  on  the  arc, 
(JFCST) 
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42-50      19        Maximum  arc  repair  cost,  i.e.,  the  cost 

of  repairing  the  arc  from  its  lower  bound 
capacity  to  its  maximum  capacity.  (RMAXCT) 

52-60      19        Maximum  arc  repair  time,  i.e.,  the  time 

required  to  repair  the  arc  from  its  lower 
bound  capacity  to  its  maximum  capacity. 
(MRTM) 

66-70      F5.3       Arc  vulnerability  parameter.   (VP) 

Last 

Card   1-2       A4        The  word  'GO'. 


Note:  Arc  data  cards  must  be  placed  in  the  following  order: 

1.  Universal  arc. 

2.  Invulnerable  and  artifical  arcs. 

3.  Vulnerable  arcs. 
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APPENDIX  B 


GLOSSARY  OF  PROGRAM  TERMS 


ARNAM(k) 

COSTMN(k) 

CRC 

CRT(k) 

CSTMNT 

DAMAGE (k) 

DAY 

FLOW(k) 

FROM(k) 

GO 

IKEEP1 

IKPOP2 

IKPOP1 

IMNCST 
IMNCUT 

INFEAS 

INVL 

IOPARC 

IRQFLO 


The  designation  of  arc  ij .  (input  data) 

Real  number  equivalent  of  MNCOST. 

Cumulative  repair  cost  for  all  interdicted  arcs. 

Cumulative  repair  time  for  an  interdicted  arc. 

Real  number  equivalent  of  MNTCST. 

Total  reduction  in  capacity  of  arc  ij  caused  by  one 
or  more  strikes  against  it. 

Counter  for  number  of  days  up  to  NUMDAY. 

Amount  of  flow  on  arc  ij . 

Location  of  beginning  node  of  arc  ij .  (input  data) 

Flag  used  to  denote  last  data  card  has  been  read. 

Used  to  save  capacity  of  arc  ij . 

Used  to  save  capacity  of  the  optimal  arc  to  strike. 

If  OPTARC  is  two-way,  saves  the  capacity  of  the 
reverse  direction  arc. 

Cost  of  the  initial  flow  in  the  network. 

Cost  of  flow  in  the  network  after  the  optimal  arc 
has  been  struck. 

A  flag  used  in  the  Out-of -Kilter  subroutine  to 
denote  infeasible  flow. 

The  number  of  invulnerable  or  artificial  arcs  in 
the  network,  (input  data) 

If  an  arc  is  two-way,  specifies  the  reverse 
direction  of  the  best  arc  to  strike. 

Integer  equivalent  of  REQFLO. 
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is  (10 

ISWTCH 
IT(k) 

IV  (k) 

JFCST(k) 

JFLOW(k) 

J2WAY 

KFLOI 

LBND (k) 

LINE 

MCAP  (k) 

MCAPU 

MCTY (k) 

MNCOST 

MNTCST 

MRTM(k) 

N 

NEWCAP (k) 

NN(k) 

NR 

NRC 

NUMARC 

NUMDAY 


Counter  for  arcs  that  have  been  struck  one  or  more 
times. 

Flag  used  to  denote  a  two-way  arc. 

Counter  for  total  strikes  on  arc  ij  during  the 
campaign. 

Counter  for  total  strikes  on  arc  ij  during  one  day 
of  the  campaign. 

Cost  of  one  unit  of  flow  on  arc  ij .   (input  data) 

Saves  flow  on  OPTARC  before  strike,  (for  output) 

Same  as  ISWTCH. 

Saves  flow  on  OPTARC  after  strike,  (for  output) 

Lower  flow  bound  on  arc  ij .   (input  data) 

Line  counter  used  in  the  output  subroutine. 

Capacity  of  arc  ij  before  next  strike. 

Saves  the  capacity  of  OPTARC  before  strike. 

Flow  capacity  of  arc  ij .   (input  data) 

Network  flow  cost  before  next  strike. 

Network  flow  cost  after  arc  ij  has  been  struck. 

Maximum  repair  time  for  arc  ij .   (input  data) 

Maximum  number  of  strikes  that  can  be  allocated 
during  any  one  day.   (input  data) 

Capacity  of  arc  ij  after  strike. 

Unique  node  number  assigned  by  Subroutine  Number. 

Counter  for  total  strikes  on  all  arcs  during  any 
one  day. 

Counter  for  total  campaign  strikes. 

Maximu  number  of  arcs  in  the  network. 

Number  of  days  in  the  campaign.   (input  data) 
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NUMNDE 
OPTARC 
PvEDU 

REPDAY(k) 

RMAXCT(k) 

RPMCST(k) 

TME(k) 

TO(k) 

TOTCST 

TOTDST 

VALUE (k) 

VP(k) 

VULCAP  (k) 

XINCRE(k) 

X2WAY 


Maximum  number  of  nodes  in  the  network. 

The  optimal  arc  to  strike. 

The  percentage  by  which  initial  flow  is  to  be 
reduced.   (input  data) 

Day  on  which  arc  ij  is  to  be  repaired. 

Maximum  repair  cost  for  arc  ij .   (input  data) 

Marginal  repair  cost  for  arc  ij . 

Marginal  repair  time  for  arc  ij . 

Location  of  ending  node  of  arc  i j .   (input  data) 

Counter  for  cumulative  campaign  strikes. 

Counter  for  cumulative  day  strikes. 

The  value  of  striking  arc  ij . 

Vulnerability  parameter  of  arc  ij .   (input  data) 

The  vulnerable  capacity  of  arc  ij . 

Marginal  reduction  in  capacity  on  arc  ij  caused  by 
an  additional  strike. 

Input  parameter  that  specifies  whether  or  not  all 
arcs  are  two-way. 
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APPENDIX  C 


PROGRAM  LISTING 


DIMENSION  XINCPE(150) 

COMMON /CATOl /I NVL ,DA Y, NUMD A Y , NUMARC , OPT ARC , L I NE t INFEAS 

C0MMON/DAT02/ST NUMNDE,MCAPUfKFLCI,NR,N 

COMMON/DATC3/!MRCtI  WNCST,  I MNCUT , X2WAY , "NCOST , MNTCST 

COMMON /DAT04/REP DAY (150) ,ARNAM ( 1 50 ) , FROM( i 50 ) ,T0 ( 150 ) 

COMMON/DAT0  5/RMAXCT(15G) ,RPMCST(15  0),PI(150) 

COMMON/CAT06/LBND(  150 )  ,NE WCP ( 1 5  0 ) , MC AP ( 15  0 ) , 1 S ( 150 ) 

COMMON/DAT  07 /MRTM( 150),NN(150),MCTY(150),JFCST(150) 

COMM0N/CAT08/FLOM 150) ,NL( 150) 

COMMON /DAT09 /VALUE (150) , PRONAM( 25 ) , JFL OW ( 150 ) 

COMMON/DAT  10 /R£QFL0,PEDU,J2WY,  ISWTCH.CRC 

COMMON /GAT  11/ I  OP ARC , TO  TOST ( 150) ,T0TCST(150) , IT (150) 

COMMON /CAT  12/ I V( 150)  , DAMAGE ( 150),VP(150) 

COMMON/DAT13/IROFLC,VULCAP(150) , CRT  (150) 

COMMON /DAT1<+/TME(  150 J  ,VULCTY(150) 

INTEGER  T0TDST,T0TCST,CAY,ARNAM,X2WAY,FP0M,T0 

INTEGER  FLOW, PI  ,OPTARC , REPDAY, S , RMAXCT ,CRC 

DATA  NAME/'GO1/ 

OPTFLO  =  0. 
REOFLO  =  0. 
IROFLO  =  0 
1  =  0 

READ  PROGRAM  TITLE  AND  PROGRAM  PARAMETERS 

READ(5,5C5)(PR0NAM(KK)  ,KK=1,20) 
READ(5,50C)N,NUMDAY,lNVLfX2WAY,REDU 
WRITE( 6,5C4) 
S  =  INVL  ♦  2 

READ  ARC  PARAMETERS 

10  1=1+1 

READ(5f51C)ARNAM( I ) ,FROM( I ) ,T0 ( I ) , LBND ( I ) , MCTY ( I) 
If  JFCSK  I  ), RMAXCT (I  )  ,  MRTMU  )  ,VP(I  ) 
IFURNAM(I)  .EQ.  NAME)  GO  TO  20 

DETERMINE  PARAMETERS  FCR  FIRST  STRIKE 

VULCAP(I)  =  MCTY(I)  -  LBND(I) 
MCAPU  )  =  KCTY(  I) 

DAMAGL(I)  =  VULCAPU  ) * ( 1 .O-EXP (-VP (  I)  )  ) 
NEWCP(I)  =  FLOAT(MCTY( I) )-OAMAGE(I  )  +  .5 
RPMCST(I)  =  (FLOAT(R^AXCT(I ) )/VULCAP(I  ))*DAMAGE(  I  ) 
TME(l)  =  (FLOAT(  MRTMd  )  ) /VJLCAP(  I)  )*CAMAGE(  I  ) 
GO  TO  10 
2C  CONTINLE 

NUMARC  =  1-1 

CALL  NUMBER  SUBROUTINE  TO  DESCRIBE  NETWORK 

CALL  NUMBER 

INITIALIZE  VARIABLES 

DO  30  K  =  1,NUMARC 

IS(K)  =  0 

IT(K)  =  0 

IV(K)  =  0 

FLOW(K)  =  0 

NL(K)  =  0 
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RfcPDAY(K)  =  0 
TOTCST(K)  =  0 
TOTDST(K)  =  0 
CRT(K)  =  0. 
VALUE(K)  =C. 
VULCTY(K)  =  0. 
XINCRE(K)  =  0. 
30  CONTINUE 

SOLVE  UNINTERDICTED  NETWORK 

CALL  0CK(NUMNDE,NUMARC,FRO",TO,JFCST,MCAP,LBND,FL0W 
1PI ,NL,  INFEAS) 
IF{  INFEAS. £0,1)  GO  TO  70 
OPTFLO  =  FLOAT(FLOWd)  ) 
RtQFLO  =  CPTFL0*(1.-R£DU) 
IROFLO  =  RnQFLO 

PRINT  PROGRAM  TITLE  AND  MODE  OF  OPERATION 

CALL  GLTPUT(l) 

CAY  =  C 

CALL  0UTPUT(6) 

FIND  MAX  FLOW, FLOW  COST, AND  PRINT  RESULTS 

CALL  FLCCST(FL0W,JFCST,NUMARC,IMNCST,INVL,1) 

CALL  0UTPUT<7) 

DAY  =  1 

NR  =  0 

NRC  =  0 

CRC  =  0 

INTERDICTION  OF  NETWORK  BEGINS  HERE 

40  NR  =  NR  ♦  1 

NRC  =  NRC  +  1 

INFEAS  =  0 

DO  45  I=S,NUMARC 

JFLOW(  I  )  =  iLOW( I) 
45  CONTINLE 

CALL  STKARC  TO  DETERMINE  BEST  ARC  TO  STRIKE 

CALL  STKARC 

IF( INFEAS. EO.l)  GO  TO  70 

OPTARC  IS  THE  BEST  ARC 

MCAPU  =  MCAP(CPTARC) 
KFLOI  =  FLCW(OPTARC) 
ISAVE  =  OPTARC 
K  =  IOPARC 

STORE  CCST  OF  REPAIR  AND  COST  OF  STRIKE 

CRC  =  CRC  +  RPMCST(OPTARC)  ♦  .5 
50  CONTINLE 

CRT(K)  =  CRT(K)  ♦  TVF(K) 
REPDAY(K)  =  DAY  +  CRT(K)  4-  .5 

COUNTER  FCR  ARCS  THAT  HAVE  BEEN  STRUCK 

IS(K)  =  IS(K)  +  1 

COUNTER  FCR  DAY  STRIKES 

IV(K)  =  IV(K)  +  1 

COUNTER  FOR  CAMPAIGN  STRIKES 

IT(K)  =  IT(K)  +  1 
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TOTDST(K)  =   IV(K) 
TOTCST(K)  =  IKK) 

C      DETERMINE  PARAMETERS  FOR  NEXT  STRIKE: 

J  =  IS(K)  ♦  1 

MCAP(K)  =  NEWCP(K) 

VULCTY(K)=MCAP( K)  -  LBND(K) 

XINCRt( K)=VULCTY(K)*<1.-EXP(-VP(K> ) ) 

NEWCP(K)=FLCAT(VCAF(K) )-XINCRE(K)  ♦  .5 

DAMAGE(K)  =  MCTY(K)  -  NEWCP(K) 

RPMCST(K)  =  (FLOAT(RMAXCT(K) J/VULC/ P(K) KXINCRE(K) 

TME(K)  =  <  FLOAT(MRTMCK)  ) / VULCAP< K ) ) »  <•  INCR£(K  ) 

C      IF  ARC  IS  2-WAY,  GO  BACK  AND  COMPUTE 
C      PARAMETERS  FOR  RETURN  ARC 

IF  ( ISWTCh.EO.l)  GO  TO  60 
K  =  ISAVE 
IShTCH  =  1 
GO  TO  5C 

60  CONTINUE 

C      PRINT  STRIKE  RESULTS 

CALL  0UTPUT(2) 

7C  CONTINUE 

IFUNFEAS  .EQ.  1)  GO  TO  130 
IF(FLOMl)  .LT.  LBND(U)  GO  TO  130 
IF(FL0W(1)  .EQ.  0)  GC  TO  140 

C  .  IF  REDL  IS  <  1, NUMBER  OF  STRIKES  IS  NOT 
C      CONSTANT.  IF  NUMBER  OF  STRIKES  =  NUMBER 
C      SPECIFIED  IN  INPUT  OR  IF  CAPACITY 
C      HAS  BEE*  REDUCED  TO  SPECIFIED  AMOUNT, 
C      DAY  IS  CONE. 

IF(REDU.LE.1.0)G0  TO  80 

IF(NR  .EQ.  N)  GO  TO  90 

GO  TC  40 
80  CONTINUE 

IF(FLOW(  D.LE.REQFLOJGO  TO  90 

IF(NR.EQ.N)GO  TO  90 

GO  TO  40 
90  CALL  0UTPUT(5) 

IF(OAY  .EQ.    NUMDAY)  GO  TO  150 

LINE  =  4  8 

C      RE-INITIALIZE  PARAMETERS  FOR  NEXT  DAY 

IOC  DAY  =  CAY  +  1 

NR  =0 

DO  110  M=S,NUMARC 

IV(M)  =  0 

TOTDST(N)  =  0 
110  CONTINUE 

C      PERFORM  REPAIR  CN  ARCS 

DO  120  M  =  S,NUMARC 
IF( IS( M)  .EQ.  0)  GC  TO  120 

MCAP(M)=FLGAT( MCAP(M) ) +( 1 . /CRT ( M ) ) *DAM AGE ( M ) 
VULCTY(M)=NCAP( M)  -  LBND(M) 
XINCRE(M)=VULCTY(M)*(i.-EXP(-VP(M) ) ) 
NEWCfM  M)=FLOAT(  MCAf(  M)  )  -  XINCPE(M)  ♦  .5 
DAMAGE(M)=MCTY(M)  -  NEWCP(M) 

RPMCST(N)  =  (FLOAT(RMAXCTIM) )/VULCAP(M) )*XINCRE(M) 
TME( M)=(FLOAT( MRTM(M) ) /VULCAPt  M )  ) *X  I  NCRE ( M ) 
CRT(M)  =  CRT(M)  -  1. 
12C  CONTINUE 
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SOLVE  NETWORK  FOR  NEXT  DAY 

CALL  OOK ( NUMNDE, NU f ARC , FROM, TO, JFCST ,MCAP, LRND,  FLOW, 
1PI  ,NL,  INFEAS) 
IF( INFEAS  .EO.  1)  GO  TC  70 

FIND   FLCk  COST 

CALL  FLCCST(FL0W,JFCST,NUMARC,IMNCST,INVL,1) 

PRINT  NETWORK  FLOW  PATTERN  AND  ITS  COST 

CALL  0LTPUT(6) 

CALL  0LTPUT(7) 

IF(FL0M  1).G7.REQFL0)GC  TO  40 

GO  TO  ICO 
13C  CALL  OLTPUTO) 

GO  TO  90 
140  CALL  OL7PUT(4) 

GO    TO    90 
150    CONTINLE 

CALL    EXIT 

500    F0RMAT(4I5,1X,F4,2) 

504  FORMAT* 1H1,10X) 

505  F0RMAT(20A4) 

510    F0RMAT(A4,2(1X,A2) ,5(1X,I9),5X,F5.3) 
END 
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c 

DEFIN 

c 

c 

NAME 

c 

c 

NODES 

c 

ARCS 

c 

I 

c 

J 

c 

COST 

c 

HI 

c 

LO 

c 

FLCW 

c 

PI 

c 

NA 

c 

INFEA 

OOK  IS  THE  FULKERSON  OUT  OF  KILTER  ALGORITHM 


SUBROUTINE  OOK ( NODES , ARC S , I , J,CCST,HI, LO, FLOW, PI 
1,NA,  INFEAS) 

DIMENSION  K 150)  ,J(150) ,COST { 15 0  ) , HI (150) , 
1L0(150),FL0W( 150) 

DIMENSION  PI (150) ,NA(150) 

INTEGER  NODES, ARCS, I ,J, C OS T, HI, LC, FLOW, PI, NA, INFEAS 

INTEGER  A, A A, N, SRC  ,SNK  ,DEL  ,  I  NF  ,C  ,  ACK  ,  CCK , Nl ,N? , 
1INC, LABEL 

DEFINITION  OF  CALLING  SEQUENCE 

USE 

NUVBER  OF  NODES 

NUMBER  OF  ARCS 

LIST  OF  FROM  (BEGINNING)  NODES 

LIST  OF  TO  (ENDING)  NOCES 

UNIT  COST  OF  FLOW  CN  ARCS 

UPPER  BOUNDS  FOR  ARCS 

LOWER  BOUNDS  FCR  ARCS 

AMOUNT  OF  FLOW  IN  ARCS 

NODE  PRICES 

SCRATCH  LIST  FCR  NCDES  LABELING 

FLAG  DENOTING  THE  CONDITION  OF  OUTPUT 

INF  =  2147483647 
AOK  =  C 

C      LOOK  FCR  AN  OUT  OF  KILTER  ARC 

AA=1 
100  Nl  =  I (AA) 
HZ    =  J(AA) 

C  =  COST(AA)  ♦  PI(Nl)  -  PKN2) 
40  IF(FLOMAA).LT.LO(  AA).OR.  (  C.LT.O.  AND.  FLOW  (  A  A  ).  LT. 
lHI(AA) ) )  GO  TO  50 

IF(FLOMAA).GT.HI(  AAJ.CR.  (  C.GT.O.  AND.  FLO  W(  A  A  )  .GT. 
1L0(AA)  )  )  GO  TO  60 
90  AA  =  AA  +•  1 

IF(AA.LE.ARCS)GO  TO  100 

C      NO  OUT  OF  KILTER  ARCS  LEFT 

INFEAS  =  0 
RETURN 

C      OUT  OF  KILTER  ARC  FCUND 

5C  SRC  =  J(AA) 

SNK  =  KAA) 

LABEL=  +AA 

GO  TO  2C0 

6C  SRC  =  KAA) 

SNK  =  J(AA) 

LABEL  =  -AA 

C      SAVE  LABELS  IF  LAST  OPERATION  WAS  INCREASING  NGDE 
C      PRICES  CN  THIS  ARC 

200  IF(AA.EO.ACK.AND.NA(SRC) .NE.O)  GO  TO  205 
DO  201  N  =1, NODES 

NA(N)  =  0 

201  CONTINUE 
AOK  =  AA 

205  COK  =  C 

NA(SRC)  =  LABEL 
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LABEL 


21C 


C 
C 


C 
C 


220 


240 


245 
250 


260 
27C 


3CC 


310 


315 
320 
340 

350 

36C 


LABEL  =  0 

DO  250  A  =1,ARCS 

Nl  =  1(A) 

IF(Nl.LT.O)  GO  TO  250 

N2  =  J(A) 

IF(NA(N1) .EQ.0.AND.NA(N2).E0.0)  GO  TO  250 

IF(NA(N1).NE.0.AND.NA(N2) .NF.O)GO  TO  245 

C  =  COST(A)  +  PKN1)  -PKN2) 

IF(NA(MI.EQ.O)  GO  TO  220 

IF(FLOW(A) .GE.HI (A) .  CR  .  (  FL  OM  A  )  .  GE  .L0( A) .AND.C.GT.O)  ) 

GO  TO  245 

NA(N2)  =  A 

GO  TO  240 

IF(FLOMA).LE.LO(A).OR.(FLOW(A).LE.HI(  A)  .AND.C.LT.O)  ) 

GO  TO  245 

NA(N1)  =  -A 

LABEL  =  1 

NODE  LABELLED,  TEST  FCR  BREAKTHRU 

IF(NA(  SMO.NE.O)  GO  TO  260 
1(A)  =  -Nl 

CONTINUE 

GO  BACK  AND  DO  MORE  LABELLING  IF  SOME  N3DE  WAS 
LABELLED  CN  LAST  PASS 

IF(LABEL.NE.O)  GO  TO  210 

RESTORE  POSITIVE  SIGNS  TO  FIRST  NODE  LIST 

DO  270  A=1.ARCS 
1(A)  =  IABS(  1(A)  ) 
CONTINUE 

IF  NO  LABELLING  DONE  ON  LAST  PASS,  GO  TO 
INCREASE  PIE 

IF(LABEL.EG.O)  GO  TO  400 

BREAKTHRU,   FIND  THE  INCREMENT 

INC  =  INF 

FOLLOW  PATH  BACK  FRCM  SOURCE 

N  =  SRC 

A  =  IABS(NA(N) ) 

IFtNA(N).LT.O)  GO  TO  315 

N2  =  I (A) 

C  =  COST(A)  -  PI (N)  +  PI (N2) 

IF(C.GT.O)INC  =  MINO(INC,LO( A)  -  FLOW(A)) 

IF(C.LE.O)INC  =  MINOUNCHI  (A)  -FLCW(A)) 

GO  TO  340 

N2  =  J (A) 

C  =  COST(A)  +  PI  (N)  -PKN2) 

IF(C.LT.O) INC  =  MINO( INCtFLGW(A)  -H I ( A ) ) 

IF(C.GE.O)INC    =    MIN0(  INC.FLOM A)    -LO(A)) 

N    =    N2 

IF(N.NE.SRC)GO    TO    310 

INCREMENTS 

A  =  IAES(NA(N) ) 

IF(NA( N).LT.O)GG  TO  360 

FLOW(A)  =  FLOW(A)  +  INC 

N  =  1(A) 

GO  TO  370 

FLOVs(A)  =  FLOW(A)  -INC 

N  =  J(A) 
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37C  IF(N.NE.SRC)GO  TO  350 

FLOW  INCREMENTED,  RETURN  TO  KILTER  TEST 

NA(N)  =  0 
GO  TO  IOC 

CHANGE  PI 

4CC  DEL  =  INF 

FIND  INCREMENT 

DO  420  A=1,APCS 
Nl  =  1(A) 
N2  =  J(A) 

IF(NA(M).EQ.0.AND.NMN2)  .EQ.OJGO  TO  420 
IF(NA(Nl).Nt.0.AND.NA(N2).NE.0)G0  TO  420 
C  =  COST(A)  +  PI  (Nl)  -  PKN2) 

IF(NA(N2)  .EQ.O.AND.FLQM  A)  .LT.HKA)  )  DEL  =MINO(DEL,C) 
IF(NA( N2).NE.0.AND.FL0M A).GT.LO(A) )  DEL  -MI NO ( DcL , -C ) 
420  CONTINUE 

IF(DEL  .NE.INF)  GO  TO  430 

IF(FLOW(AA).EQ.LOl AA ) . CR.FLOW( AA ) . EQ.H I (  AA ) ) GO  TO  425 

INFEASIBLE  SOLUTION 

INFEAS  =  1 
RETURN 

INCREASE  PI 


425  DEL  =  lABS(COK) 
430  DO  450  N  =ltNODES 

IF(NA(N).EQ.O)  PUN) 
45C  CONTINUE 


=  PI(N)  +  DEL 


GO  BACK  TO  KILTER  TEST 

GO  TO  100 

END 
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SUBROUTINE  FLOCST{ FLCW ,COST , NUMARC »MINCT , I NVL , JUN I V ) 
DIMENSIGN  FLCW( 150) ,C0ST(150) 
INTEGER  FLCW,COST 

MINCT  =  0 
IUNV  =  INVL 

IF  JUNIV  =  Ot  USE  CCST  CF  UNIVERSAL  ARC 

IFUUNIV  .EO.  0)  IUNV  =  1 
DO  10  I  =  IUNV,NU^ARC 
10  MINGT  =  MINCT  +    FLOW (I)*COST ( I ) 

RETURN 
END 


SUBROUTINE  NUMBER 

COMMON/DAT01/INVLtOAYfNUMDAY,NUMARCf OPTARCtLINE, INFEAS 

CGMMON/DAT02/S,NUMNQE,MCAPU»KFL0I,NR,N 

COMMON/ CAT Ct/REPDAY( 150) ,ARNAM<150) , FROM ( 1 50 ) t TO ( 150 ) 

COMMON/DAT07/MRTN<150) , NN( 150) , MCT Y ( 150 ) , J FCST < 150 ) 

INTEGER  TOTDSTtTOTCSTtDAYf ARNANf X2WAYfFR0Mf TO 

NUMBER  THE  NODES 

NUMNDE  =  0 

DO  30  L=ltNUMARC 

NN(NUMNCE  +  1)  =FRCM(L) 

FROM(L)    =    NODNUM(FROML)  ,NN) 

NUMNDE  =MAXO<NUMND£,FRCP(L) ) 
3C  CONTINUE 

DO  40  L  =1,NUMARC 

NN1NUMNDE  ♦  1)  =  TC(L) 

TO(L)  =  NODNUM(TO(L) ,NN) 

NUMNDE  =  MAXO(NL^NDEfTO{L) ) 
4C  CONTINUE 

CHECK  FCR  DEAD  NODES 

DO  48  I  =1,NUMNDE 

DO  42  L  =1,NUMARC 

IFtFROP  (L)  .EO.DGO  TO  44 
42  CONTINUE 

WRITE(6,522)NN<  I) 
44  CONTINUE 

DO  46  L  =1,NUMARC 

IF(TG(  D.EC.1IG0  TO  48 
46  CONTINUE 

KRITE(6,530)NN(  I) 
48  CONTINUE 

RETURN 

522  FGRMATdHl,  22HNG  ARC  BEGINS  AT  NODE  ,  A4/1HI) 
530  FURMATdHl,  22HN0  ARC  ENDS  AT  NODE    ,  A4/1H1) 
END 


FUNCTION  NCDNUM( II  I ,NN) 

DIMENSION  NN(1) 

COMMON/ 0 AT 02/S, NUMNDE, MCAPU,KF LCI ,NR,N 

NCDNUM  =  NUMNDE  ♦  1 
IFCNUMNCE  .EO.ORETURN 
CO  10  I  =  i, NUMNDE 
IF( I  II  .EQ.NM I ) )  NCDNUM  =  I 
1C  CONTINLE 
RETURN 
END 
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SUBROUTINE  STKARC 

COMMCN/DAT01/INVL,DAY,NUMDAY,NUMARC,CPTARC,LINf ,  INFEAS 
COMMON/CAT02/S,NUMNDE,MCAPU,KFLCI,NR,N 
COMMON / C AT03/NR C, I NNC ST, I MNCUT,X2 WAY, MNCOST, MNTCST 
COMMON /DA T04/REP DA Y( 150) ,ARNAM( 150 ) , FROM ( 1 50 ) , TO ( 1 50 ) 
C0MMGN/CAT05/RMAXCT( 150),RPVCST(150),PI(150) 
COMMON /DA T06/LBND( 1 50)  ,NE WCP( 1 5  0 ) , *CAP < 150 ) , I S ( 1 50 ) 
C0MM0N/DAT07/MPTM( 150) ,NN( 150) ,  MCTYQ50) , JFCST(150) 
COMMGN/DATC6/FL0K( 150) ,NL( 150) 

COMMON/CAT09/V41UE(150) ,PR0NAM(25) t JFL0W(150) 
C0MM0N/DAT10/RE0FLn,REDUtJ2WY,ISWTCHfGPC 
CGMMON/DAT11/IOPARC,TOTDST( 150 ) , TOTCST ( 1 50 ) , I T ( 15  0 ) 
COMMON /DAT  12/1 V( 150)  ,DAMAGE(150),VP(150) 
COMM0N/DAT13/IR0FLG,VULGAP(150) ,  CRT (150) 
CCMMGN/DAT14/TMEU50)  ,  VULCTYU50) 
INTEGER  TOTDST f TOTCST, DA Y,ARNAM,X2 WAY, FROM, TO 
INTEGER  FLOW, PI , CPTARC , REP  DAY, S , RMAXCT , CRC 

C      INITIALIZE  PARAMETERS 

MNTCST  =  C 

MNCOST  =  0 

CSTMNT=0. 

COSTMN=C. 

UPVAL  ■  -1CCOO000. 

C      SOLVE  FCR  LAST  FLOW  PATTERN 

CALL  OCMNUMNDEt  NUNARC,FRGMfTO,JFCST,MCAPf  LEND,  FLOW, 
1PI ,NL, INFEAS) 
IF(  1NFEAS,EG.1)RETURN 

CALL  FLCCST(FLOW. JFCST , NUMARC , MNCOST , I NVL , 0 ) 
DO  100  I  =  S, NUMARC 
IF(FLOWU  ).EQ.O)GQ  TO  100 
IKEEP1=MCAP(I ) 
MCAP( I )=NEWCP( I  I 

C      SOLVE  INTERDICTED  FLCW  PATTERN  FOR 
C      SELECTED  ARCS  AND  FIND  FLOW  COST 

CALL  OCK ( NUMNDE, NUMARC, FROM, TO, JFCST, MCAP  ,  LBND,  FLCW , 
1PI,NL,  INFEAS) 
IF(  INFEAS. EQ.DRETURN 
CALL  F LCCST( FLOW, JFCST, NUMARC, MNTCST, I NVL, 0) 

C      FIND  VALUE  OF  STRIKING  ARCS  AND 
C      SAVE  MAXIMUM  VALUE 

MCAPU  )  =  IK£EP1 

CSTMNT=MNTCST 

COSTNN  ONCOST 

VALUE (  I  )=(CSTMNT-COSTMN)*TME(I  )+RPMCST(I) 

IF<  VALUE(  I  J.LT.UPVADGC    TO    90 

UPVAL=VALUE( I ) 

C  OPTARC    IS    THE    BEST    ARC    TC    STRIKE 

CPTARC=I 

90    CONTINUE 

100    CONTINUE 

DO  52  I  =  S, NUMARC 

IF (FROM (OPTARC) . EO . T0(  I ) . A  NO . TO ( OPTARC ) . EO. FROM ( I ) 
l.AND.ARNAM(OPTARC) .EO.ARNAMd ) )GC  TO  54 

52  CONTINUE 

53  CONTINUE 
J2WY  =  1 

C      10PARC  IS  THE  REVERSE  DIRECTION  ARC 

IOPARC  =  CPTARC 
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GO  TO  56 
54  J2WY  =  2 

IOPARC  =  I 

ISWTCH  =  J2hY 

IF(X2WAY.tQ.l)  GO  TO  53 
56  CONTINLE 

IKPOP2  =  KCAP( IOPARC) 

IKPOP1  =  NCAP(CPTARC) 

MCAP(OPTARC)  =  NEWCPICPTARC) 

KCAP(ICPARC)  =  NEWCP( ICPARC) 

STRIKE  BEST  ARC  AND  FIND  FLCW  CCST 

CALL  OCMNUMNDEf  NUPARCfFROKt  TOf  JFCSTtMCAPt  LBNCf  FLCW, 
1PI,NL,INFEAS) 
IF(  INFEAS.EO.DRETURN 

CALL  FL CCST (FLOW, JFCST, NUMARC , IMNCUT , I NVL ,  1) 
MCAP(CPTARC)  =  IKPCP1 
PCAP(ICPARC)  =  IKPCP2 
RETURN 
END 
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SUBROUTINE  OUTPUT* J) 

COMMON/DAT  01 /I NVLtD AY, NUMDAYtNUMARCt OPT  ARC tLTNEt  TNFEAS 
COMMON/DAT02/S,NUMNQE,MCAPU,KPLOI ,NR,N 
C0MMON/DATO2/NRCIMNCST. IMNCUT , X2WAY, MNCOST, MNTC ST 
C0MMON/DATOA/RFPOAY(150l ,ARNAM(150) TFROM( 150  )  ,T0( 150) 
COMMON/ CATC 5/RMAXCT( 150) ,RPMCST( 150) , PI ( 150) 
C0MMCN/DATO6/LBNDU50)  »  NEWCP  ( 1  50  ) ,  MC  AP  (  1  50  ) ,  IS  (150) 
COMMON /O AT 07/MRTM( 150) ,NN(15C) ,MCTY(150) , JFCST(15C) 
COMMON/ DATC8/FL0W( 150) ,  NL( 150) 

COMMON/DAT  0  9 /VALUE (150) , PRONAM(25)  » J  FLOW (150) 
C0MM0N/DAT10/REQFL0,REDU,J2WY, ISWTCH,CRC 
C0MM0N/DAT11/I0PARC,T0TDST(150 ) , TOTC ST( 1 50 ) , I T( 1 5C ) 
COMMON/DAT  1 2/ I V( 150), DAMAGE (150) ,VP(150) 
C0MMON/DAT13/IRQFLOtVULCAP(150) , CRT (150) 
INTEGER  TOTCSTtTOTCSTtDAYt ARNAM,X2WAY, FROM, TO 
INTEGER  FLOW, PI , OP TARC , RE PDA Y , S , RMAXCT,CRC 

IFUoGTol  )  GO  TO   20 

WRITE(6,506) (PRONAM( I) ,1=1 ,20) 

IF(REDUoLTolo  )G0    TO    10 

WRITE(6,507)NUMDAY,N 

LINE    =    8 

RETURN 
10  WRITE(6,510)N,NUMDAY,FL0W(1) ,  I  ROFLO,  I ROFLO,N 

WRITE(6,504) 

LINE  =  0 

RETURN 
20    IF(LINEoNEoA8)    GO    TO       30 

WRITE(6,5C4) 

LINE    =    0 
30    IF(LINEoNEoC)GO    TO       40 
40    GO   T0(100, 100, 150, 160, 210, 50, 180), J 
50    IARC    =    0 

IF(DAYoEOoO  ) WRITE (6, 50 2) 

IF(DAYoEQo0)G0    TO    70 

WRITE(6,519)DAY 
60    LINE    =    0 
70    WRITE(6,508) 
80    IARC    =    IARC    +    1 

IFUARCoGToNUMARC)     GO    TO    90 

IF(FLOW(  IAROoEQoOo  ANDoDAYoNEoO)    GO    TO    80 

LINE    =    LINE    +    1 

II    =    FROM( IARC) 

IJ    =    T0(  IARC) 

WRITE(6,50  9)ARNAM( IARC),NN( II) ,NN( IJ) ,LBND(IARC) , 
lMCAP(IARC)  , JFCST(IARC) ,RMAXCT( I  ARC ) , MRTM( I  ARC ) , 
2FL0WUARC)  ,VP(  IARC) 

IF(LINEoNEoA9)    GO    TO    80 

WRITE(6,504) 

GO  TO  60 
90  CONTINUE 

LINE  =  0 

RETURN 
100    CONTINUE 

IF(LINE    +60GE0A8)    GO    TO    200 

IF(LINEoNEoO)    GO    TO    120 
130    CONTINUE 

WRITE(6,5C?)DAY 

WRITE(6,50O) 

LINE    =    LINE    +    A 
120    CONTINUE 

K    =    OPTA&r 

III=FROM(K) 

JJJ=TO(K) 

WRITE (6, 501  )NR,NN( I  II ) ,NN( J J J)  , MCAPU , MC AP( K)  , 
UFLOW(K)  ,KFLOI,REPDAY(K)  ,CRC,FLOW(  1)  ,  IMNCUT 
2,IS(K) 
130    LINE    =    LINE    +    1 

RETURN 
150    WRITE(6,513)DAY 
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LINE  =LINE  +  8 

RETURN 
160  WRITF(6,512)DAY 

LINE  =  LINE  +  8 

RETURN 
180    WRITE<6,515)     IMNCST ,FL0W( 1 ) 

IF(DAYoE0o0)G0    TO    1^0 

IF(FLOW(l)oGToREOFLO)GO    TO    190 

WRITE(6t 516JDAY 
190    CONTINUE 

WRITE(6,504) 

LINE    =    0 

RETURN 
200    LINE    =48 

GO   TO      20 
210    CONTINUE 

WRITE(6,520)DAY 

WRITE(6,521)DAY 

DO    220    JARC=S,NUMARC 

KK    =    FROM(JARC) 

KJ    =    TO(JARC) 

IF(TOTDST( JARC)oEQoOoANDoTOTCST(JARC )oEQo0) 
1    GO    TO    220 

IF(FLOW(JARC)oEOoO)GO    TO    220 

WRITE(6,52  2)ARNAM( J  ARC) t NN ( KK ) , NN( KJ ) ,T0TDST ( J ARC ) i 
1TOTCSTC JARC) 
220    CONTINUE 

WRITE(6,523)NR,NRC 

RETURN 


500  F0RMAT(T30, • ARC 


CAPACITY' ,T48, f ARC 
i 


FLOW 


,T6l, 
i 


•TO  BE' 


1T71, 'CUMo • ,T77,«  MAXIMUM'  ,T 87,' NETWORK' ,T96 , ' NUMBER • , / 
2,T10,'MISSICN' ,T28, 'BEFORE  AFTER' , 4X, ' BEFOR E  AFTER' 
3,2X,'REPAIRED',2X,'REPAIR',1X, 'NETWORK' , 4X , ' FLOW • , 4X , 
4«0F  ARC  ,/,T10,  'NUMBER',  2X,  'FROM  TO' , 2X , « STRI KE ' , 2X , 
5'STRIKE' ,3X, 'STRIKE  STRIKE • , 2X, 'ON  DAY ', 4X, 'COST ', 3X, 
6 "FLOW ,6X, 'COST' ,4X ,« STRIKES • ) 


501  FORMAT (lOXt I3,5X,A2,3X,A2, 2X,  2 ( I  5, 3X ) , 1 X, 2 ( I  5 , 3X ) , 
114, 5X, 15, IX, 17, IX, 19, 3X, 14) 

502  FORMAT ( «0' ,19X,' INITIAL  FLOW  PATTERN  FOR  LOC  NETWORK' 
It/) 

503  FORMATCO' ,8X, 'STRIKES  TARGETED  FOR  DAY', 13, 'o') 

504  F0RMAT(1H1,10X) 

506    FORMATdH    ,  28  X,20A4  ,///  ) 


507    FORMATdH 
112H    DAYS, 


,32X,     23HCAMPAIGN    WILL       LAST    FOR  ,13 

WITH       ,13,'     MISSIONS    PER    DAYo  '  ) 


508    F0RMAT(44X, 'LOWER' ,6X, 'UPPER' , 13X, «MAX< 


5X, 'MAXo • , 


It  •ARC ,/,20X, 'ARC ,4X, 'FROM' ,5X,'TQ',6X, «FL0W',7X, 


15X 


2' FLOW' , 5X,  'FLOW'  ,4X,« RE  PAIR'  ,3 X, 'REPAIR'  ,4X, 'ARC  ,7X, 
3'VULNERo  «,/,20X,  'NAME' ,3X, 'NODE' ,4X, 'NODE' ,5X, 'BOUND' , 
46X,« BOUND' ,4X, «C0ST',5X, 'COST' ,5X,'TIME',5X, 'FLOW' , 
56X, 'PARAMo ' ) 


509    F0RMAT(2H       ,18X, A4 ,4X, A2 , 6X , A2 , IX, I  8, 5X, I  8 , 1 X, 17, 
12(4X,I5) ,1X,I8,6X,F5o3) 


510    F0RMAT136X,  'M/^XIMUM    NUMBER    OF    STRIKES     EACH 
1, 13, /,36X, 'NUMBER    OF    DAYS    IN    THE    CAMPAIGN 


2//,36X, 'NOTE 
3/,36X,'0NE    OF 


4« 


J3C 
II 


UNINTEPDICTED 
t 


MISSIONS    WILL 
THE    FOLLOWING 


FLOW    ( 
5«REDUCED   TO    ',14,'    UNITSo' 
6'BEGINNING    CF    EACH',/,40X, 


DAY  :  « 

:'  ,13, 
CONTINUE  EACH  DAY  UNTIL' 
OCCURS  :',/,36X, 


7«    UNITSo',  /,36X,  M3I 
8,'TARGETEDo' ) 


13 


14, •    UNITS)     IS' ,/,40X 


,/,36X, • (2)     FLOW    AT    THE     • 
•DAY     IS    REDUCED    TO    ' 
•     STRIKES    HAVE    BEEN 


,14, 
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512  FORMAT (1H  , ///10X, 27HFLOW  TOTALLY  STOPPED  ON  DAY   , 
1110,///) 

513  FORMAT (1H  , ///10X, 30HMIN IMUM  THRU-PUT  CANNOT  BE  MET, 
1"  ON  DAY' t 110,///) 

515  FORMAT( «0» ,//,20X, 'FLOW  COST   :',I8,/,20X, 
1« THROUGHPUT  : • , 18) 

516  FORMAT( «0'  , //,20X, 'NO  STRIKES  REQUIRED  ON  DAY*, 13, 'o') 

519  FORMAT* «0' ,19X,' ALL  ARC  REPAIR  SCHEDULED  FOR  DAY', 13, 
1'  HAS  BEEN  COMPLETEDo  ',/, 20X, 'RESULTING  FLOW  PATTERN' 
2, •  IS  AS  FOLLOWS: • ,/) 

520  FORMAT( ' 0'  ,/,20X, 'RECAP  FOR  DAY' , I  3 ,':«,/, 20X, 
I • ***************** i  ) 

521  FORMAT! '0' ,45X, 'TOTAL  NUMBER     TOTAL  NUMBER' ,/ ,21 X , 
l1  ARC ,21X,  'ARC  STRIKES      ARC  STRI KES « , / , 2 2X, • NAME ' , 
24X, 'FROM' ,5X, «T0«,5X, 'FOR  DAY' , I  3, 6X , • FOR  CAMPAIGN') 

522  F0RMAT(2H   , 20X, A4, 5X, A2 , 6X , A2 , 9X, 14, 12X , I  4) 

523  FORMAT{ '0' , /,22X, • NUMBER  STRIKES  THIS  DAY       :', 
119, /,22X, 'NUMBER  STRIKES  THIS  CAMPAIGN  :',I9) 

END 


47 


APPENDIX  D 


COMPUTER  OUTPUT  FOR  SAMPLE  PROBLEM 
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